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1. Introduction 



QCD dynamics at low momenta can be described in terms of a chiral effective theory, which is 
formulated by assuming chiral symmetry and its spontaneous breaking. The effective Lagrangian 
is expanded as [[l]-^ 

where 

= ^Tr [d^U^d^^U] - |Tr [e'^^^^'^U + t/l-^^g-'^/^f] , (1.2) 
= I'^-^/- (1-3) 

i 

U S SU(A'f) parametrises the pseudo-Goldstone bosons degrees of freedom; £ and F are the (in- 
finite volume) quark condensate and pseudoscalar decay constant in the chiral limit ' ; ^ is the 
A'f X Nf quark mass matrix. For completeness, the dependence on the vacuum angle 6 is also 
shown. At NLO, 10 terms appear for N{ = 2, whereas ^j^'^^ contains 12 terms for Nf = 3, with 
corresponding couplings 

{%} ^ li=LJ,hi=i.3 [Nf = 2], (1.4) 
{-^■l ^ L,=i..io,//,=i..2 [A^f = 3]. (1.5) 

A more comprehensive introduction on Chiral Lagrangians has been presented by J. Bijnens at this 
conference [Q]. F, £, {/,} ({L,}) are the so-called Low-Energy couplings (LECs): they parametrise 
the low-energy dynamics, which is not determined by symmetries and can in principle be derived 
from QCD by means of a non-perturbative method. Once the LECs aie known, chiral perturbation 
theory becomes a powerful predictive framework for investigating low-energy properties. Lattice 
QCD represents the ideal tool to match QCD with the chiral effective theory and to extract LECs in 
a reliable way; unquenched simulations are now reaching masses and volumes where such a match- 
ing can be performed, and many efforts have been made with the final goal to keep all systematic 
errors under control. 

On a finite volume V = L^T with L » 1 /Aqcd, different chiral regimes can be distinguished. 
Approaching the chiral limit by keeping Mj^L » 1 defines the so-called /7-regime, where the power 
counting in terms of the momentum p and the quark mass m is given by 

mr^p^, \/L,\/Tr^p. (1.6) 



In this regime, the chiral effective theory in a finite box looks very much like in infinite volume: 
finite- volume effects are exponentially suppressed by factors ~ exp(— M^jL), while mass-effects 
are dominant. 

Alternatively, one can approach the chiral limit by keeping /i = niLV < 0(1); this defines 
the so-called £-regime |^, where the Compton wavelength associated to the pseudo-Goldstone 

'in the following F and L will refer to the case Nf — 2; for A'f = 3 the notation Fq, Zq will be used. 
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bosons is larger than the linear extend L of the box, Mj^L < 1 ^. In this case the power counting is 
reorganised such that 

m~£^ 1/L, l/r~e. (1.7) 

It follows that mass effects are suppressed, while volume effects are enhanced and become polyno- 
mial in L^^. One of the consequences of the rearrangement is that, at a given order in the effective 
theory, less LECs appear with respect to the /j-regime: the NLO predictions are less contaminated 
by higher order effects, making the £-regime particularly advantageous and clean to compute the 
leading order couplings F and Z. Furthermore, in the £-regime topology plays a relevant role [^. 
Observables can be defined at fixed value of the topological charge, and the dependence on this 
charge should also be well reproduced by the effective theory. Therefore topology provides a new 
variable in this regime, in addition to the mass and the volume. In the infinite volume limit, this 
dependence is expected to vanish. 

The possibility of matching the chiral effective theory with QCD in different regimes and by 
means of different observables opens the chance for several independent determinations of LECs 
-with different higher order corrections and systematic errors - from which a precise and first- 
principles determination of the LECs can be achieved in the near future. 

2. Lattice determinations of LECs in the p-regime 

2.1 Nf = 2 

In QCD with two flavours, the chiral expansion of the pion mass and pion decay constant at NLO 
in the isospin limit is given by |^ 



M2 . /M^ 

with 



M- = (2.3) 
The scale-independent LECs It, , Z4 are defined as 

^3,4=ln(Ai,4/M^)M=139.6MeV (2.4) 

The values of these parameters extracted from phenomenological analysis are [|2[, |8|- 10 1 



F = 86.2(5) MeV, /3 = 2.9(2.4), Z4 = 4.4(2). (2.5) 
The full NNLO expression for those quantities is also known (see [§]). 

In many recent lattice studies with A'^f = 2 dynamical fermions the quark mass dependence of Mj^ 
and Fji has been studied in order to extract F, £, 1-^, I4. The main characteristics of the different 
computations are summarised in table |l], while the results are listed in table ^; notice that the quark 



^The case Mj^L < 1 and T ^ L, corresponding to the so-called 5-regime, will not be considered in this discussion. 
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condensate is given in the MS scheme at energy scale 2 GeV. 

The first attempt to estimate directly 73,4 from unquenched lattice calculations has been performed 
in [jlTI], where two-flavour QCD with Wilson quarks (with and without non-perturbative on-shell 



0{a) improvement) has been simulated using DD-HMC algorithm [12, 13]. The authors pointed 
out that the extraction of 74 is problematic, since the obtained results are not stable by adding higher 
order terms in the NLO fit of the pseudoscalar decay constant. The error quoted in table § for 73 is 
statistical only. 

ETM collaboration is carrying on an extensive project with Nf = 2 Wilson twisted mass quarks 
at maximal twist. First results with light quarks have been presented in [p4||. Updated results have 



been presented at this conference 1 15, 16] and are reported in table g; the data have been corrected 
for finite-size effects using the extended LUscher formula |p^]. The quoted errors are statistical, 
0{a^) effects and finite-size effects respectively. Notice that in this study the scale has been fixed 
through the physical values of Mj^ and F„. See |[l5|, [l^ for the complete set of ensembles produced 
by the collaboration. 

QCDSF and UKQCD are performing Nf = 2 simulations using Wilson fermions with non- 
perturbative on-shell 0{a) improvement. Several lattice spacings and volumes have been simu- 
lated (see [18] for the full list). The data ai^e con^ected for 0{a^) effects and for finite-size effects 
estimated from [17]. The scale in this case has been determined from the axial coupling gA- 

Finally, JLQCD presented at this conference the first results obtained with overlap Dirac op- 
erator and Iwasaki gauge action [19, 2C]. Two kinds of finite volume corrections have been taken 
into account: the usual finite volume effect computed according to [17], and the correction due to 
fixed topology, according to [21 1. Large differences have been observed by adding 2-loop terms 
in the chiral fit, with respect to a NLO fit; two different NNLO chiral fits have been tried, and the 
relative difference considered as systematic uncertainty (second error in table §). Here the scale 
has been fixed by the condition ro = 0.49 fm [^]. 

All these results fairly agree among themselves; although systematic errors need to be better 
quantified, we expect them to be under control in the next years, so that the path towards a reliable 
computation of LECs of two-flavour chiral theory looks very promising. 

The phenomenological impact of the lattice computations of 73 , 74 has been discussed recently 
by H. Leutwyler [^ ; in particular he considered the 5-wave pion scattering lengths (/ = 0,2), 
which can be related to 73,74 by means of low-energy theorems [^. 



2.2 Aff = 2 + 1 

Lattice results obtained with three dynamical flavours can be matched with Nf = 3 chiral perturba- 
tion theory in order to extract Lo,Fo,{Li}. In particular, in the three-flavour chiral effective theory, 
the quark-mass dependence of pseudoscalar decay constants at NLO contains the LECs L4 and L5, 
while pseudoscalar masses are sensitive to the combinations (2L8 — L5) and {2L(, — L4) More- 
over, the same LECs can be extracted from a partially quenched setup [Q |^] (see also [p6|]). 
The current status of phenomenological determinations of the L, has been summarised by Bijnens 
in this conference [Q] ; from a NLO fit one infers 

L4 = 0; Lfi^O; 10^ -15 = 1.46; 10^ -18 = 1.00. (2.6) 
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Authors 


Dirac op. 


gauge act. 


a (fm) 


r/a, L/a 


^;r,min 

(MeV) 


Del Debbio etal|ll] 


Wilson 

Wilson 0{a) impr. 


Wilson 
Wilson 


0.0717(15) 

0.0521(7) 

0.0784(10) 


64, 32 
48,24 


403 
381 
377 


ETM [111, 16] 


Wilson Twisted Mass 


Sym. tree 


0.0858(5) 
0.0858(5) 
0.0657(11) 
0.0657(11) 


48,24 
64, 32 
48,24 
64, 32 


300 
300 
420 
300 


QCDSF/UKQCD [|l|] 


Wilson 0{a) impr. 


Wilson 


-0.08 
-0.08 
~0.07 


48,24 
64, 32 
64, 32 


450 
340 
440 


JLQCD [19, 20] 


Overlap 


Iwasaki 


~0.12 


32, 16 


290 



Table 1: Parameters of recent A^f = 2 simulations in the p-regime with light quarks. The second column 
refers to the Dirac operator, the third to the gauge action adopted. In the sixth column, M;^ min indicates the 
lightest pion mass reached in the simulations. 



Authors 


F 

(MeV) 


ri^S(2 GeV) 
(MeV)^ 


h 


U 


Del Debbio etal [11] 






3.0(5) 




ETM [T|, 16] 


85.98(7)(21)(35) 


266(6)(0)(6) 


3.44(8)(26)(6) 


4.61(4)(3)(7) 


QCDSF/UKQCD [18] 


79(5) 


273(12) 


3.49(12) 


4.69(14) 


JLQCD [[H, H] 


78(3)(1) 


242(6) (6) 


2.9(4)(2.6) 


4.3(5)(2) 



Table 2: Summary of LECs obtained from Nf = 2 simulations. 



The parameters related to recent N{ = 2 + \ lattice simulations are summarised in the top part of 
table H; the corresponding results for the renormalised {2L^ —L5), {2L(, — L4), L4, L5 evaluated at 
the scale Mp = 770 MeV are collected in the bottom part. 

The MILC collaboration presented the first results obtained with staggered fermions and fourth 
root prescription in 2004 ||27|]; several updates followed (see for instance the contribution at the 



2006 lattice conference [28]) and new preliminary results were presented at this conference 



The chiral fits are performed using (partially quenched) rooted staggered chiral perturbation [ |30| ] 
including analytic NNLO and NNNLO terms. By setting the scale through Fji they obtain for the 
LO constants 

= 1.052(2) (+^), £p(2GeV) = [278(l)(+^)(5)MeV]3, (2.7) 
F^/Fo = \.2l{5){t'!), Sr(2GeV) = [242(9) (4) MeV]^, 

where Fq,'Lq refer to the 3 flavour chiral limit, while F2, £2 are obtained by sending mu^md —>■ and 
Ms —>■ phys. The quoted errors are statistical and lattice-systematic; for the condensates there is 
an additional uncertainty coming from perturbative renormalisation. 
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Authors 


Dirac op. 


gauge action 


a (fm) 


L(fm) 


MpsMn (MeV) 


MILC[|^ 


impr. staggered 


Sym. 1 loop 


0.06-0.15 


2.4-3.4 


200 


RBC-UKQCDI33] 


Domain Wall 


Iwasaki 


0.11 


1.8-2.6 


330 


PACS-CS[f5|l 


Wilson 0{a) impr. 


Iwasaki 


0.09 


2.9 


210 



Authors 


(2L8-L5)-103 


(2L6-L4)-10^^ 


u ■ 10^ 


L5 • 10^ 


MILCI29] 


0.3(1)(1) 


0-3(1) (^3) 


0.1(3)(+^) 


1.4(2) (+^) 


RBC-UKQCD[||] 


0.247(45) 


-0.002(42) 


0.136(80) 


0.862(99) 


PACS-CS[35] 


-0.23(5) 


0.10(4) 


-0.02(11) 


1.47(13) 



Table 3: On the top: parameters of recent Nf = 2 + \ simulations in the /j-regime. On the bottom: results 
for the LECs evaluated at the scale Mp = 110 MeV. 



RBC and UKQCD collaborations presented this year first results for light meson masses and 
pseudoscalar decay constants computed with Domain Wall Dirac operator and Iwasaki gauge action 
pT]]; in this first work, with L ~ 1.8 fm and pseudoscalar masses up to 400 MeV, they observed 
that a simultaneous NLO fit of pseudoscalar meson masses and decay constants fails. New results 
have been presented at this conference, with a larger volume and lighter quark masses [32|, 33] (in 
particular, Mps^mm — 330 MeV for the full case and M/>5 min — 250 MeV in the partially quenched 
setup). By observing that three-flavour fits are not safely applicable in the mass range up to the 
kaon mass, they also performed partially quenched SU(2) x SU(2) fits, where only two quarks are 
treated as light. To compare the two fits, they derived from both cases the constants 1-^,1 a using 
NLO matching between two- and three-flavour LECs [^, finding compatible results within the 
errors. 

Finally, PACS-CS collaboration is undertaking Nf = 2-\-\ simulations with 0{a) non- pertur- 
batively improved Wilson quarks and Iwasaki gauge action, adopting DD-HMC algorithm; prelim- 
inary results were presented at this conference |34, 35 1. 

Given the poor knowledge of systematic errors, a precise comparison of the lattice results at 
this point is maybe premature. 



3. Determinations of LO LECs in the e-regime 

Lattice simulations in the £-regime require small quark masses, and in general the preservation 
of chiral symmetry at finite lattice spacing would be highly desirable: the spectral gap of the 
Dirac operator is then bounded from below and this ensures stability of the dynamical simulations. 
However it is a matter of fact that Ginsparg-Wilson fermions are still very demanding from the 
numerical point of view, since the construction of the corresponding Dirac operator is very costly 
with respect to other regularisations. 

On the other hand, for the Wilson Dirac operator the spectral gap is not guaranteed to have a lower 



bound and recent studies established an empirical stability condition [36| 



a 

m > mmin; rUi-nin °^ —f=. (3.1) 



6 



Determining QCD Low-Energy Couplings from lattice simulations 



Silvia Necco 



In order to reach the region mZV < 1, still unrealistically large (with the present computing re- 
sources) lattice extents would be needed ^ . 

The spectral gap of the twisted mass Dirac operator is also bounded from below: an attempt to 
investigate the £-regime with this discretisation has been presented at this conference [38|. 

Another reason that makes the e-regime particularly challenging is that in the region where 
miy ~ 0(1) one expects large fluctuations associated to low-modes wavef unctions, which can 
induce large uncertainties on the observables. These can be however substantially reduced by 
adopting the so-called low -mode averaging technique [|^, 40 1. 

At leading order in the £-expansion, the partition function of the effective theory reads 



SU(N) 



dUo exp 



Tr(e'*^/^^[/o + e-'*^/^[/,^^' 



(3.2) 



where Uq represents the collective constant field associated to the zero modes. The partition func- 
tion at fixed topology v is obtained by Fourier-transforming in Q : 



(!¥' 



U(N) 



J[/o(det[/o)''exp 



^Tr (^^Uq + [/J^^ 



(3.3) 



3.1 Quark condensate from finite-size scaling 

We consider for simplicity a quark mass matrix proportional to the identity, ^ 
and introduce the dimensionless variable }i 
% 

>: H 

Em 



ml/^p A'^f > 2, 

ml}/; the quark condensate for /i ^ 1 behaves like 



r(At) = — ^ In iT^ 



(3.4) 



As expected, vanishes in the chiral limit since spontaneous symmetry breaking does not occur 
at finite volume. Analogously, at fixed topology one has p] 



with Xv 



2{Nf+v 



■At + . 



(3.5) 
(3.6) 



where the infrared divergence proportional to 1/m is due to zero-modes contribution. From these 
expressions it is clear that, even if spontaneous symmetry breaking does not occur in a finite vol- 
ume, the formation of a quark condensate £ leaves signs in a finite box and the corresponding 
finite-size scaling is predicted by the chiral effective theory. 

In particular, is exphcitly known up to NLO. The O(e^) contributions are entirely given by 
one-loop corrections to the chiral condensate 

t2 



1 + 



(3.7) 



where jSi is a known universal shape coefficient [41 and L = V^^^. The important point is that 
NLO corrections still contain only LO LECs, a typical fact which occur in the £-regime, as already 



■'For 0(a)-improved Wilson fermions there are first indications that the width of the spectral gap distribution has a 
different scaling Win, although definitive conclusions are still missing. 
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pointed out in the introduction. 

From the QCD side, given a lattice regularisation which satisfies the Ginsparg-Wilson relation, the 
quark condensate at fixed topology can be expressed as 



pv(A) 



2m / dl- 



Vm 



(3.8) 



where Pv(A) is the spectral density of the Dirac operator and l/m divergence is exactly the one 
founded also in the effective theory. Once renormalised with the logarithmic-divergent renormal- 



isation constant Z5, Xv can be matched with Xv of eq. 3.6. Notice that at finite mass Xv contains 
also additive ultraviolet divergences which must be removed in order to perform a finite-size scal- 
ing study. Being topology-independent, a convenient solution for this purpose is to consider the 
combinations {Xvi — Zv2)> which have an unambiguous continuum limit at finite quark masses. 
Several quenched studies have been performed in order to extract £ from a finite-size scaling study 
^ [0-0- The left side of fig. |l| shows the quantity {X2 — Xs) computed in [£7| for two symmetric 
volumes V = {1.5 fmY and V = (2.0 fm)^, as a function of (mV). The fact that {X2 - Xs) /{mV) 
does not depend on the volume within the statistical errors indicates that LO scaling is verified, 
and there is no sensitivity to NLO corrections. These data have been obtained using low-mode 
averaging technique. 

The chiral effective theory gives non-trivial parameter-free predictions in the chiral limit, also 
known as Leutwyler-Smilga sum rules ; for instance from eq. |3|^ one obtains 



Xv,{^) - Xviil^) 



^1=0 



(Vl-V2)(V3+M-)(V4+M-) 
(V3-V4)(Vi+M-)(V2+M-)' 



(3.9) 



A comparison with the lattice data p7| ] in the quenched approximation is shown in the right part 
of fig. H, for several values of the topological indices and for three different lattices. The quark 
mass dependence of the ratios is weak, and the data refer to a very light quark mass corresponding 
to /I ~ 0.07. The good agreement supports the fact that the topology dependence of those ratios is 
well reproduced by the (quenched) chiral effective theory. 



3.2 Meson correlators 

The pseudoscalar decay constant can be extracted by matching current correlators in the £-regime 
with the NLO predictions. For example, for the left-handed cunent correlator at fixed topology 



(3.10) 



one obtains the following prediction at NLO 

72 



C"{xo) = Tr 

^Xo\ 1 



1 

12 



(3.11) 
(3.12) 



*In the quenched hmit Zgff (K) diverges logarithmically with L and hence can be defined only at a finite volume. 
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0.2 
_ 0,15 

T?. 0.1 
>< 

0.05 


3x10-" 



a = 0.125 fm. L— 1.5 fm 
a = 0.125 fm, L=2.0 fm 



100 



200 
mV 



300 



400 



1.4 



1.2 - 



theory 



0.8 



0.6 - 



0.4 



(1-3) 



s 

(1-4) 
(1-3) 



(1-5) 
(1-3) 



, a — 0.125 fm, L=1,b fm 
, a = 0,125 fm. L=2,0 fm 
, a = 0,09 fm, L=l,5 fm 



Figure 1: On the left: finite-size scaling of differences {X2 — X^)^ computed at two different symmetric 
volumes V = (1.5 fm)'* (black empty squares) and V = (2.0 fm)^ (red filled points) p?]]. The data refer to 
the quenched approximation, with overlap Dirac operator. On the right: the ratios {Xv, — XvT)/{Xvi ^ Xv^) 
for three different lattices considered in p7|] compared with theoretical predictions of eq. 3.9 



T" are the SU(A'^f) generators, Zj is the renormalisation constant associated to the left current, j3 



and /:()() depend on the geometry and 'Lv{pL) is defined in eq. Also in this case, up to NLO 
only the LECs £ and F appear. Current correlators in the £-regime have been investigated in many 
quenched studies ^ using Ginsparg-Wilson fermions. In [40] the JlJi correlator has 



been computed both in the e- and /j-regime: here it has been shown that the quenched value of F 
obtained in the /7-regime after a chiral extrapolation is in agreement with the one extracted from 
the £-regime. 

JLQCD presented at this conference preliminary unquenched results obtained with overlap 



Dirac operator and Iwasaki gauge action, A^f = 2 []55|]. They computed PP, SS, VoVo and AqAq 
correlators at a lattice spacing a = 0.11 fm, V = 16^ x 32 and fixed topology v = 0. Fig. ^ 
shows the time-dependence of PP and AqAq correlators, for a/Wsea = fl^'^vai = 0.002 corresponding 
to /I = 0.556. By simultaneously fitting pseudoscalar and axial correlators they obtain 

F = 87.3(5.5)MeV, £^(2 GeV) = [239.8 (4.0)MeV]^ (3.13) 

Pseudoscalar and current correlators have been recently computed in the effective theory at NLO 
also for non-degenerate quark masses in the full and partially quenched scenarios, in particular for 
the case where all quarks are in the p- or in the £ -regime and for the mixed case, where mvaiSV < 1 



and mseaSV' » 1 [Q^, |57p. These results may be very useful for the future matches of lattice QCD 
with the chiral effective theory. 

An alternative method to extract pseudoscalar decay constant from two-point correlators in the e- 
regime has been proposed in [^]: by matching residuals of 1/m^ poles with the expectations of 
the chiral effective theory, one can extract F by means of zero-modes correlation functions at fixed 
non-zero topology. 
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0.01 



0.008 



lattice ma=0.002 
ChPT (E=[227MeV]^ F =89.5MeV.) 



0.002 




10 



15 



20 



25 



30 



0.0001 



8e-05 



6e-05 - 



with LMA ma=0.002 
without^LMA )77n=0.002 
ChPT (Z=[260MeV] , f_=89.5MeV.) 



4e-05 - 



2e-05 - 



10 



15 



20 



25 



30 



Figure 2: Pseudoscalar and axial vector correlators computed by JLQCD with overlap Dirac operator, Nf = 2 
3.3 LO LECs from eigenvalue distribution 

At LO in the £-expansion, the partition function at fixed topology is equivalent to the one of a chiral 



Random Matrix Theory (RMT) []59|-|62|]; it follows that RMT reproduces the same microscopic 
spectral density p5v(C)M) the chiral effective theory, in terms of two dimensionless variables 
= Xiy and n = mLV, where A are the eigenvalues of the Dirac operator. Moreover, it is possible 



to extract the probability distributions of single eigenvalues [63-65]: 



k=l 

By matching the low-lying spectrum with these expectations 



(A«FVrL2) 



(3.14) 



(3.15) 



one can then extract the low-energy constant £. 

The low-lying spectrum of the Dirac operator in the quenched approximation has been matched 
with the predictions of RMT in many works [^-|6^. In [ ]68| ] a detailed comparison has been 
made, finding a good agreement for L > 1.5 fm. Starting from this knowledge, several unquenched 



computations have been performed in the past years [70-73] : their relevant parameters are sum- 
marised in table ^ JLQCD/TWQCD [ ]72| ] performed two separated sets of simulations, in the e- 



and the /7-regime; in the last case they studied the sea quark mass dependence of £ and verified that 
its extrapolation to the chiral limit is compatible with the result obtained in the £-regime. Moreover, 
the sea quark mass dependence of ratios of eigenvalues for A/f = 2, v = is such that close to the 
massless limit reproduces the Nf = 0,v = 2 results (flavor-topology duality), while at heavy masses 
becomes compatible with Nf = 0, V = 0. Notice that in [ |70| , [tT] ] quark masses are still in a regime 
where mLV is considerably larger than 1 ; nevertheless the observed quark-mass dependence of the 
extracted low-energy constant is negligible within the statistical errors. In figure ^ the results of 
JLQCD/TWQCD are shown as an example. Ratios of eigenvalues {l^ic)/{Q) are compared with the 
expectations of RMT; in particular the central plot corresponds to A^f = 2, v = 0,/i = 0.556(16); 
the left plot refers to A/f = 0, V = and the right to A/f = 0, V = 2. The physical results for the quark 
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Authors 


Dirac op. 




a (fm) 


L/a,T /a 


mLV 


V 


1^(2 GeV) 
(MeV)3 




overlap 


2 


0.15 


10,10 


-2-5 


0,1 


282(10) 




chirally impr. 


2 


0.11-0.13 


12,24 


-5-15 


0,1 






overlap 


2 


0.11 
0.11-0.12 


16,32 
16,32 


0.556(16) 
-7-30 




0,2,4 


251(7) 


1^3,74] 


fixed point 


2+1 


0.13 


12,12 


- 1.5 (u,d), - 13 (s) 


0,1,2 





Table 4: Summary of simulation parameters of recent Nf = 2 simulations performed to extract the quark 
condensate by matching low-lying spectrum with RMT predictions. 
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© 

"®" 


-*- 
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^4/2^ 


-er -sr 


-e- -e- 




3/2 4/3 
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-e- 



Figure 3: Ratios of eigenvalues (Q) / (^/) for several k,l computed in [ [72| ] compared with expectations of 
RMT. 



condensate in the MS scheme at 2 GeV are reported in the last column of table 
An important issue concerning these computations is the estimation of higher order corrections: 
the matching with RMT is valid only at leading order and hence no control of higher order effects 
is possible. Assuming that finite-volume corrections for £ are the same as in the chiral effective 
theory, one would obtain a systematic error of 31(11) MeV respectively on r'^^(2GeV)'/^ for 
[f70|]([72]). Another open question is that, while it is clear how the spectral density is renormalised 



[]36|], this is not the case for the individual eigenvalues. 

This framework can be extended such that the spectrum of the Dirac operator is sensitive also to 
the pseudoscalar decay constant F at LO in the chiral effective theory; a proposal has been made 
in [^-^]. One considers the eigenvalue equations of the Dirac operator in presence of external 
constant Abelian gauge potentials with couplings /^iso.i 7^ Miso.2: 

D,,2wia = [^(A) + /Atiso,i,27o] wii = i^i'Mi- (3-16) 

Translated into the chiral effective theory language, it leads to a modified partition function, which 
at LO in the £ -regime reads: 



'V 



dUoidetUoYexp 

U(N) 



(3.17) 



For |[7lL p3p the computation of the renormalisation constant Zs is still missing. 
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with 

B = diag(Atiso,ilA',,Miso,2liV2); Nf = Ni+N2. (3.18) 
Note that the imaginary nature of the chemical potential guarantees a positive definite fermion 



determinant. Moreover, generalised Leutwyler-Smilga spectral sum rules can be derived [ |79| ] 



The same partition function can be obtained starting from a chiral Random Two-Matrix Theory [ |78[ ] 



(R2MT), from which it is possible to compute the spectral correlation functions in the microscopic 
limit Pj^^''^j)({Ci}ni {?2}/t) for every Nf = N\ +N2 and («,/:). The predictions from R2MT -or 
equivalently from the LO chiral Lagrangian- are now formulated in terms of the rescaled variables 
(^1 = X{LV , (^2 = ^i^y and }JLf = mfiy. Moreover, the coupling with the isospin chemical potential 
introduces a new dependence on the combinations jUiso,i,2^V^, opening the possibility to extract 
F by matching appropriate correlation functions measured on the lattice with the expectations of 
the LO effective theory. A typical example is the mixed correlation function 

plf;liX,M = - - A^'"')) - - - (3.19) 

n m n m 

The effect of having a non-zero chemical potential can be observed for instance in fig. ^ [[77||, where 
the prediction of the microscopic correlation function is plotted in the dynamical Nf = 2 case, with 
Ni = N2 = I, Miso,! = — Miso,2 = Miso- The variable (^2 = 4 is fixed, and the i^i -dependence is plotted: 
for /liso = the curve would show a delta-function in correspondence of l^i = 4; for /^iso 7^ 
one observes a non-zero width. The quark masses in this case would change the height of the 
curve. As for the ordinary random matrix theory, it is theoretically possible to extract individual 
eigenvalue probabihty distributions in terms of all density correlators. Truncated expansions have 
been investigated and presented at this conference [^]. 

First numerical tests of the method have been performed in jTSj |7^ ], using unimproved staggered 
quarks at coarse lattice spacing, both in the quenched and the dynamical Nf = 2 case. 



A pilot study with overlap fermions, Nf = 2 [ ]81p , has been presented at this conference | |82| ] ; a lattice 
spacing a ~ 0.13 fm has been adopted, with L/a = T /a = \2 and mLV ~ 4. In this work, the case 
A'^i = Ni = 2, N2 = 0, /Xiso,i = 0, /Xiso,2 = Miso has been considered: it corresponds to a partially 
quenched situation, with two sea quarks at chemical potential /liso.i = and two valence quarks 
coupled to /^iso. The fit of the integrated correlator yields the results lM^{2 GeV) = [234(4) MeV]^ 
(with finite- volume corrections estimated from chiral effective theory) and F = 101(6) MeV (with 
no finite-volume correction). 

This method to extract F can be relatively cheap compared for instance to the computation of 
current correlators; however, as already mentioned, NLO are not under control and -in absence of 
theoretical developments in this sense- simulations at several volumes would be needed in order to 
estimate finite-size effects. 



4. Conclusions 



Lattice QCD simulations are now reaching quark masses and volumes where the matching with 
the chiral effective theory can be performed in a reliable way, keeping the different sources of sys- 
tematic errors under control. Many results are already available in the /7-regime for Nf = 2; from 
the NLO quark-mass dependence of Mji and Fji it is possible to extract the LLCs Zs, 74, F and £. 
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Figure 4: The microscopic correlation function p'',1^'-',)(^i, ^2) in eq- 3.19 with /Xiso.i = — Miso.2 = M 



5,(1,1)^ 



m„ ^/EV =5, = 4 in the figure) , as a function of ^1 ((^+) [Q. The red dashed curve represents the 
case (with the delta-function omitted), while the black solid curve corresponds to j^l^iJ^^V ~ 0.05. 



Mis 



Many efforts have been spent to evaluate discretisation errors and finite-volume corrections. Cur- 
rent simulations are now reaching pion masses of 300 MeV; in order to control uncertainties from 
higher order corrections in the chiral effective theory, lighter quark masses would be needed. With 
the recent algorithmic improvements and increasing computational resources, this is a realistic goal 
for the next years. 

New results are also available from A^f = 2 + 1 simulations: besides the MILC results, RBC/UKQCD 
and PACS-CS presented preliminary studies at this conference, obtained with domain wall and 
0{a) improved Wilson fermions respectively. From the quark-mass dependence of pseudoscalar 
masses and decay constant it will be possible to give reliable estimations of Fo,Lo, L4,L5,L^,L^. 

In the e-regime, £ and F can be extracted from observables which at NLO are not contami- 
nated by higher order LECs. One can obtain fully independent determinations with respect to the 
ones achieved in the p-regime; their agreement is non-trivial test of QCD in the chiral regime. New 
unquenched results have been presented at this conference, both for the low-lying Dirac spectrum 
and meson correlation functions; although the £-regime simulations are challenging from the com- 
putational point of view, they provide an alternative and complementary approach to the standard 
p-regime determinations. 

In this review only QCD LECs have been considered; however, lattice computations can be 
used for other constants which encode low-energy properties, for instance to evaluate couplings of 
the chiral weak Hamiltonian or of the Heavy Meson chiral effective theory. 
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